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The optimally doped and underdoped region of the t — J model at large N (N is the number 
of spin components) is governed by the competition of d-wave superconductivity (SC) and a d 
Charge-Density Wave (d-CDW). The partial destruction of the Fermi surface by the d-CDW and 
the resulting density of states are discussed. Furthermore, c-axis conductances for incoherent and 
coherent tunneling are calculated, considering both an isotropic and an anisotropic in-plane momen- 
tum dependence of the hopping matrix element between the planes. The influence of self-energy 
effects on the conductances is also considered using a model where the electrons interact with a 
dispersionless, low-lying branch of bosons. 

We show that available tunneling spectra from break-junctions are best explained by assuming 
that they result from incoherent tunneling with a strongly anisotropic hopping matrix element of 
the form suggested by band structure calculations. The conductance spectra are then characterized 
by one single peak which evolves continuously from the superconducting to the d-CDW state with 
decreasing doping. The intrinsic c-axis tunneling spectra are, on the other hand, best explained by 
coherent tunneling. Calculated spectra show at low temperatures two peaks due to SC and d-CDW. 
With increasing temperature the BCS-like peak moves to zero voltage and vanishes at T c , exactly 
as in experiment. Our results thus can explain why break junction and intrinsic tunneling spectra 
are different from each other. Moreover, they support a scenario of two competing order parameters 
in the underdoped region of high-T c superconductors. 



I. INTRODUCTION 



One important topic in high-T c superconductors is the 
question of how many order parameters are needed for 
a proper description of the optimally doped and under- 
doped cases. One scenario assumes that only the super- 
conducting order parammeter is relevant. The decrease 
in the transition temperature T c is then caused by fluctu- 
ations of its phase and the pseudogap is locally just the 
superconducting gap. A second scenario assumes that 
the physics in the underdoped and optimally doped re- 
gion is mainly determined by the competition of the su- 
perconducting order parameter with a second one in the 
particle-hole channel. Examples could be the antiferro- 
magnetic, s- and d- charge density wave or stripe order 
parameters. 

Many experiments such as angle-resolved photoemis- 
sion [1] or tunneling in break-junctions [2-5] suggests that 
there is only one energy scale related with the gap and 
that this scale increases monotonically with decreasing 
doping. Recent intrinisc c-axis tunneling spectra in sev- 
eral cuprates [6-13] seem to modify this picture. Opti- 
mally doped or underdoped samples show at low temper- 
atures two peaks for positive or negative voltages. The 
peak at larger voltages stays essentially at the same posi- 
tion, but becomes broader with increasing temperature. 



With decreasing doping it moves towards larger volt- 
ages. Though this peak behaves similar to the one seen 
in tunneling in break-junctions, it has recently been ar- 
gued that heating effects could seriously affect this peak 
[14,15]. The peak at smaller voltages moves towards zero 
voltage with increasing temperature, looses hereby in- 
tensity and vanishes at T c . Heating effects should be 
unimportant for the behavior of this peak. Intrinisc tun- 
nel spectra of this kind have been found both for double 
layer and as well as single layer materials [12]. On the 
other hand, strongly overdoped samples show only one 
sharp peak with properties as expected from BCS theory 
[6]. 

It is tempting to associate the two peaks observed in 
the optimally doped and the underdoped region with the 
SC and the pseudogap, as has been done in some of the 
above references. In the following we will investigate 
whether the widely accepted t — J model supports such a 
picture. To this end we will present calculations for the 
conductance based on a t-J model where the two spin 
components have been generalized to N components and 
the leading diagrams at large N are taken into account. 
As discussed in detail in Ref. [16] the phase diagram in 
this limit has a quantum critical point (QCP) separating 
at T=0 the normal phase at large dopings from a d-CDW 
state at lower dopings if superconductivity is omitted. 



Allowing also for superconductivity the QCP separates a 
pure superconducting state from a ground state contain- 
ing both superconductivity and a d-CDW. The proper- 
ties at optimal doping and in the underdoped regime are 
mainly determined by the competition between supercon- 
ductivity and the d-CDW. This model thus represents an 
example for the above second scenario. 



II. DENSITY OF STATES AND FERMI SURFACE 
IN THE PRESENCE OF SC AND D-CDW 



The d-CDW order parameter, appropriate for the t-J 
model at large N, is given by $(k) = — i/ZNc^ J(k — 
q)(cJ lcr Cq+Q ( j). J is the Heisenberg coupling, cJ, c are cre- 
ation and annihilation operators for electrons under the 
constraint that double occupancies of lattice sites are ex- 
cluded, N c is the number of primitive cells, (...) denotes 
an expectation value, and Q is the wave vector of the 
d-CDW. Keeping only the instantaneous term in the ef- 
fective interaction, the superconducting order parameter 
is A(k) = l/2iV c £ (J(k - q) - V c (k - q))(5 qT c_ qi ). 
As shown in Ref. [16] it is in general necessary to 
include the Coulomb potential Vc in order to stabi- 
lize the d-CDW with respect to phase separation. In 
the presence of the two order parameters the operators 
(c^ *, c_k,j,, +, c_k-Q,j.) are coupled leading to the 
following Green's function matrix [16] 
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e(k) is the one-particle energy, e(k) = — (St + 
aJ)(cos(k x ) + cos(ky)) — 2t'Scos(k x )cos(k y ) — /i, with 
a = ~L/N c J2q C0S (qx)f(c(q)). f is the Fermi function, S 
the doping away fom half-filling, /i a renormalized chem- 
ical potential, t and i! are nearest and second-nearest 
neighbor hopping amplitudes, z a complex frequency, and 
k = k Q. 

Expressing the expectation values in the order pa- 
rameters by Go and using Eq.(l) one obtains two cou- 
pled equations for the order parameters. In the in- 
teresting doping region the order parameters have d- 
wave symmetry, $(k) = $7(k),A(k) = A7(k), with 
7(k) = (cos(k x ) — cos(k y ))/2. The equations for A and 
$ read then 
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k ct=l 



( 7 2 (k)(A 2 + $ 2 ) - (E a (k) + e(k))(E a (k) + e(k)) . (3) 



Vc is the expansion coefficient of the Coulomb poten- 
tial in the d-wave channel with the basis function 7, and 
E a (k) are the four poles of Gq(z, k) in the z-plane. At 
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FIG. 1. Order parameters $ and A as a function of tem- 
perature in units of t for the dopings 8 — 0.114 and 0.078. 

zero temperature $ decreases monotonically with in- 
creasing S, whereas A first increases, passes then through 
a maximum at So and finally decreases again, as shown 
in Fig.l of Ref. [17]. Fig.l shows the temperature de- 
pendence of $ and A, calculated with t'/t = —0.35 and 
J/t = 0.3. The energy unit is t. A repulsive nearest- 
neighbor Coulomb interaction was also included with 
Vc/t — 0.06. According to Fig.l the temperature de- 
pendence of the order parameters is mean-field like suffi- 
ciently away from So- Near So the two order parameters 
strongly interact with each other. For instance, for the 
slightly underdoped case of So — 0.114, the increase of A 
at low temperatures is accompanied by a decrease of $ 
so that the "total gap" VA 2 + <i> 2 is rather constant at 
low temperatures. 

Fig. 2 contains quasi-particle densities for S = 0.114 
and T = 0. The thin dotted line denotes the density 
for vanishing order parameters. It shows a logarithmic 
divergence due to the van Hove singularity. The latter 
lies for the chosen parameters just below the Fermi en- 
ergy which corresponds to zero energy. The long-dashed 
curve decribes the case where the correct finite value for 
$ has been used, but A has been put to zero. This den- 
sity shows a strongly asymmetric gap structure with a 
strong peak on the low and a weaker, splitted peak at 
the high-energy side. A closer analysis shows that the 
lower peak of this dublett comes from k-points near the 
antinodal points X and Y. States near these points are 
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FIG. 2. Density of states for $ = A = (thin dotted line), 
$ / 0, A = (long-dashed line), $ = 0, A / (dash-dotted 
line), $ / 0, A / (solid line) for T = and 5 = 0.114. 

coupled by the d-CDW and their energies are shifted by 
the formation of the d-wave gap. This explains why this 
peak moves upwards (see the dashed and solid lines in 
Fig. 2) in the presence of an additional BCS gap. In 
contrast to that, the upper peak of the dublett originates 
from k-states on the boundaries of the reduced Brioulin 
zone near the points (jr/2 — S,tt/2 — 8) and equivalent 
points with 8 << ir/2. Their energies are mainly deter- 
mined by the one-particle energies near this point relative 
to the Fermi energy and thus less sensitive to the forma- 
tion of the gap. This explains why the position of this 
higher peak of the dublett is unchanged by the BCS gap, 
see the dashed and solid lines in Fig. 2. The density is 
everywhere nonzero in the d-CDW state, in particular, 
at the Fermi energy. The asymmetry of the density with 
respect to zero energy is caused by the asymmetric bare 
density due to the proximity of the van Hove singularity. 
The finite density of states at the Fermi energy allows to 
lower further the ground state energy by introducing a 
superconducting gap. The density of states becomes then 
strictly zero at the Fermi level and the additional d-wave 
gap is rather symmetric with respect ot the Fermi energy. 
The resulting density of states for this case $ ^ 0, A ^ 
is given by the solid line in Fig. 2. Comparing this line 
with the dot-dashed line, which corresponds to a pure SC 
state $ = 0, A ^ 0, one recognizes that the inner part of 
the gap structure looks like a reduced SC gap with weakly 
developed edges, at least on the low-energy side. Fig. 3 
shows the density of states for 8 = 0.114 and three dif- 
ferent temperatures. These temperatures are low enough 
so that the main gap edges do not change much because 
V <E> 2 + A 2 is nearly constant. However, the opening of 
the SC gap in the inner part of the gap structure can 



FIG. 3. Density of states for T = 0. (solid line), T = 0.007 
(dashed line), T = 0.011 (dotted line) and S = 0.114. 



clearly be seen as a function of temperature. 
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FIG. 4. Density of states for three different dopings at 
T = 0. 



Figs. 2 and 3 illustrate the fact that the total gap is 
not just one single d-wave gap with an amplitude given 
by the square root of the sum of the square of the two 
gaps. Instead, the SC and CDW gaps interact with each 
other, however in such a way that their individual struc- 
tures can still be seen in the density of states. This is 
also apparent in the density plots for three different dop- 
ings in Fig. 4. The upper and lower panels illustrate the 
difference between the density of states for a SC and a d 
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CDW gap, respectivley. Some features of the individual 
gaps are still present in the middle panel of Fig. 4 which 
describes the case of coexisting SC and d CDW gaps. We 
find in contrast to Ref. [18] that our self-consistently cal- 
culated order parameters yield for all considered dopings 
and hopping parameters densities where the SC gap lies 
inside the d-CDW gap. 

One important feature of the coexising SC and d-CDW 
state is that the two gaps have different locations in k- 
space: The CDW gap mainly resides near the X and Y 
points, and the SC gap near the diagonal. This becomes 
clear by looking at the Fermi lines as a function of <I> in 
the absence of superconductivity. Since a finite $ implies 
a doubling of the elementary cell we have plotted Fermi 
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FIG. 5. Fermi lines at T = for four values for $ in the 
absence of superconductivity. 



lines in Fig. 5 in the reduced zone scheme, e.g., the new 
Brillouin zone (BZ) is bounded by a straight line between 
the points X = (ir, 0) and Y — (0, 7r). The black line cor- 
responds to $ = and describes the usual normal state 
Fermi line in the reduced Brillouin zone for the parame- 
ters t'/t = 0.35 and J/t = 0.3. For $ = 0.02 the Fermi 
line consists of a long arc around the diagonal ending at 
the new boundary of the BZ and two pieces near the X 
and Y points. This means that the region around the 
hot spots becomes first gapped. Increasing $ only the 
arc around the diagonal survives and becomes shorter. 
The finite density of states at the Fermi energy in the 
CDW state is due to this arc. Allowing also for SC the 
arc becomes gapped and the Fermi line shrinks to one 
point on the diagonal. The coexistence of SC and CDW 
thus becomes possible because the CDW state can lower 
the free energy by introducing a SC gap along the arc. 



III. CONDUCTANCE 

Using lowest-order perturbation theory in the inter- 
layer hopping the quasi-particle c-axis current J between 
superconducting layers (SIS junction) is given by [19] 



J(V) = 




f(u + eV))A n (k% u + eV)A 11 (qq', u). (4) 

In Eq.(4) Tk q denotes the hopping matrix element be- 
tween states with momenta k and q in adjacent layers 
and <C ... S> an average over quenched disorder. / is the 
Fermi function, V the applied voltage, and, using the 
Nambu representation, An the spectral function of the 
element 11 of the 2x2 Green's function matrix. The mo- 
menta in the above formula refer to the original (large) 
BZ. The differential conductance G(V), which is the main 
quantity of interest in the following, is defined as the first 
derivative of J with respect to V. 

We make the following Ansatz for the averaged squared 
tunneling matrix element 

« Tk q T k *, q , »= ti 7 (k)7(q) 7 (k')7(q'Wk- q ,k'- q ' 

(*5k, q + ff(k-q)). (5) 

The form factors 7(k) determine which electrons in the 
BZ are mainly involved in the tunneling process. Results 
from band structure calculations [20] suggest the Ansatz 
[21] 

7 (k) = 1 — u + u/2\cosk x — cosky\, (6) 

with < u < 1. The parameter u interpolates between 
the isotropic case u = and the strongly anisotropic case 
u = 1. The latter is typical for tunneling within a double 
layer of CuOi planes, whereas tunneling between layers 
in different elementary cells may include also an isotropic 
component. The first term in the brackets in Eq.(5) ac- 
counts for coherent scattering with strength a. The sec- 
ond term in the brackets describes incoherent scattering 
where g is a smooth function of the momentum. One may 
distinguish two cases for the momentum dependence of 
g. In the case of strong, localized scatterer g may as- 
sumed to be completely independent on momentum. If 
long-ranged random fields are present g is large (small) 
mainly at small (large) momentum transfers. The mo- 
mentum dependence of g thus can be modelled by 

g(k)=g-expHk\ 2 /Ai), (7) 

where the momentum A interpolates from isotropic to 
forward scattering, described by large and small values 
for A, respectively. 

Since the case of coherent scattering can be obtained 
from that of incoherent scattering by replacing g by 
47T/A 2 and taking the limit A — > we will first consider 
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incoherent scattering. Using Eq.(5) and restricting the 
momenta to the reduced BZ because of the cell doubling 
due to the d-CDW we obtain for the tunnel current 

P t 2 . _ r°° 
J ( y ) = r^E^ 2 ( k )^) / <M/M - /( w + eV )) 

^ c ~ J — oo 

kq 





Ai+2i,i+2i' (k, w + eV")Ai + 2i',i+2i(q,u;) 

it' =0,1 

+ p(k - q - Q)(An (k, w + eV)A 33 (q, w) 

+ i43s(k,a; + eV)i4ii(q,w))). (8) 

The tilde on the momenta indicates that these momenta 
lie in the reduced BZ. The spectral functions Aij can be 
assumed to be periodic with resepct to the reduced BZ. 
The same is true for the form factors 7 but not for the 
function g which originates from impurity potentials. To 
make the expression for J(V) independent of the choice 
for the reduced BZ we also translate back the momentum 
appearing in g to the reduced BZ. Eq.(8) then becomes 

p t 2 . _ r°° 

J{ y) = -^E^ 2 ( k ^ 2 ® / <m/m - /( w + eV » 

kq 

g(k - q) ((An (k, uj + eV) + A 33 (k, w + eV))(4n (q, w) 
+ A 33 (q, w)) + 2A 13 (k, w + eV)A 31 (q, w)) . (9) 

The first contribution in the big parantheses describes 
s-wave, the second one d-wave scattering. Their relative 
importance is controlled by the parameter p c in the func- 
tion g. If g is independent of the transferred momentum 
Eq.(8) simplifies to 

j(y) = ^f du(f(u) - f(u; + eV))p(uj + eV)p(uj), 
n J-00 

(10) 

with the weighted density 

~P^) ^ ^ET 2 (k)(^n(k^) + A 33 (k^)). (11) 



c - 

k 



For a SIN junction one usually assumes that g is in- 
dependent of momentum. Its current is then obtained 
from Eq.(10) by identifying one of the two densities with 
that of the normal metal pu which can be assumed to be 
constant. The resulting conductance of a SIN junction 
becomes then at not too high temperatures 



G SIN (V) = e -^p{eV). 



(12) 



Since we will mainly consider SIS junctions in the follow- 
ing conductance will always refer to SIS junctions unless 
it is stated otherwise. 



The simplest case of incoherent scattering corresponds 
to u = and A = 00, i.e., where the averaged tunneling 
matrix element is independent of all momenta. Using the 
densities of Fig. 4 the resulting conductance curves are 
shown in Fig. 6. In the pure superconducting state (up- 
per panel) the conductance shows a broad peak near the 
gap 2A which decays rapidly towards larger but rather 
slow towards smaller energies. The conductance is pos- 
itive for all frequencies, especially also above 2A, which 
is intimately connected to the presence of the large and 
rather constant density of states outside of the gap re- 
gion. In the d-CDW case (lower panel in Fig. 6) the con- 
ductance curve has two peaks. The higher and dominant 
one is due to the CDW gap 2<F In contrast to the su- 
perconductor the d-CDW state has (neglecting the tiny 
BCS-gap at the doping 6 = 0.077) a finite density of 
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FIG. 6. Incoherent c-axis conductance calculated for u=0, 
A = 00, T=0, J=0.3, t'=-0.35, and three different dopings. 

states at and near the Fermi energy along the arcs. The 
folding of these states with one of the CDW edges causes 
the lower peak at about the energy <!>. Well above 2$ the 
curve is again rather constant and partly slightly neg- 
ative. In the weakly underdoped case (middle panel in 
Fig. 6) the conductance curve shows essentially two peaks. 
They arise due to the folding of the large CDW shoul- 
der at negative frequencies with the BCS peak and the 
two splitted CDW shoulders at positive frequencies, re- 
spectively. The BCS gap itself is seen only as a broad 
and weak structure at low energies. Some of the features 
in Fig. 6 agree with the tunneling experiments, e.g., the 
monotonic increase of the dominant high-frequency peak 
with decreasing doping and the appearance of more than 
one peak in the optimally and underdoped cases. How- 
ever, several details of these curves are not found in the 
experiments: The peak in the overdoped case is much too 
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broad compared to that in the intrinisic c-axis tunneling 
spectra of Ref . [6] , the lower peak in the underdoped cases 
is caused by the relaxation of electronic states around the 
CDW gap to states near the arcs or the nearby BCS gap 
and thus does not approach zero at T c as in intrinsic tun- 
neling spectra. Similar conclusions are reached, follow- 
ing Eq.(12), by comparing the experimental conductance 
curves of SIN junctions [3,4] with the densities of Fig. 4. 
The monotonic increase of the distance between the two 
main peaks with decreasing doping occurs in both cases 
but the experiment does not show the asymetry of the 
theoretical conductance curve in the underdoped case as 
well as the additional structures obtained in the region 
of coexisting SC and d CDW. 

Things change substantially if the form factor 7(k) 
with a non-zero value for u is taken into account. Assum- 
ing still a momentum-independent function g the tunnel 
current is now to be calculated from the weighted density 
p as given by Eqs.(lO) and (11). Fig. 7 shows p{w) for 
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FIG. 7. Weighted density p as a function of frequency for 
T=0, J=0.3, t) = 0.004, t'=-0.35, and three different dopings 
8. 

the extreme anisotropic case u = 1 for three different 
dopings. Most of the background contribution to the 
density has been removed. In the over- and under-doped 
cases (upper and lower panels in Fig. 7) p consists of just 
two rather symmetric peaks with respect to uj = which 
are related to the superconducting and d-CDW gaps, re- 
spectively. In the slightly underdoped regime, where SC 
and d-CDW coexist, p is still rather symmetric with re- 
spect to uj — and consists of four peaks. The peaks at 
large energies are dominant, their frequencies are roughly 
given by ±2y/$ 2 + A 2 , i.e., they describe the "total" gap 
of the two components. The two weaker and less pro- 



nounced peaks at smaller frequencies are related to the 
SC gap, which can be concluded from their temperature 
dependence and magnitude of their energies. The low 
intensities of these peaks can also be easily understood: 
The BCS gap resides on the arcs near the diagonal. The 
form factor 7 with u = 1 suppresses heavily the tunneling 
of states in this region. It is also interesting to note that 
only the lower peak of the splitted high-energy d-CDW 
edge in Fig. 4 for 6 = 0.114 survives in the corresponding 
weighted density in Fig. 7. This can easily be understood 
by noting that the lower (higher) peak of the dublett is 
due to k-states near the antinodal (nodal) points and 
thus unaffected (suppressed) by the form factor. 

The curves for p look in many respects similar to the 
experimental SIN conductance curves. In both cases, the 
spectra are dominated by two pronounced peaks lying 
rather symmetrically with respect to u> — and whose 
separation increases monotonically with decreasing dop- 
ing. These peaks evolve in p very smoothly from a SC 
to a d CDW state passing through a region where both 
order parameters coexist. The agreement can be further 
improved if one introduces a phcnomcnological damping 
in the theoretical curves. The peaks are then broadened 
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FIG. 8. Incoherent c-axis conductance calculated for u=l, 
A = 00, T=0, t'=-0.35, and three different dopings. 

and the low-energy structures for S = 0.114 become in- 
visible. The dip on the high-energy sides of the main 
peaks in the experimental spectra is, however, missing in 
p indicating the presence of self-energy effects beyond a 
constant damping. 

Performing the frequency integral in Eq.(lO) with the 
weighted densities p one obtains the curves of Fig. 8. The 
main effect of the inclusion of the anisotropic form fac- 
tors 7 with u = 1 is the suppression of the small quasi- 



6 



particle excitations near the nodal regions. This means 
in the overdoped case 8 = 0.178 that the slowly decay- 
ing tail of the main peak towards lower voltages seen in 
the upper panel of Fig. 6 is substantially suppressed mak- 
ing the peak much sharper. The dip above 2A is also 
more pronounced than in Fig. 6 and the conductance as- 
sumes (small) negative values over a wide region towards 
higher voltages. The reason for this negative resistance 
becomes clear from a comparison of Figs. 4 and 7. Most 
of the rather constant background density in Fig. 4 has 
been removed by the anisotropic form factor. However, 
just this constant background density is responsible for a 
positive and structureless conductance outside of the gap 
region. Similar considerations apply to the underdoped 
case (5 = 0.077. The anisotropic form factor sharpens 
up somewhat the high voltage peak and suppresses the 
lower peak at around $ because the states near the arcs 
can no longer contribute much. At the same time the 
conductance shows a well-pronounced dip above 2$ with 
large negative values due to the eliminated background 
density of states. Similar statements hold for the slightly 
underdoped case. Here the lower peak at around A + $, 
which was in Fig. 6 still the strongest one, is suppressed 
but still visible. 

The above calculations show that an incoherent tun- 
neling model with a momentum independent function g 
is not able to produce a peak in the coexistence region 
which moves towards zero voltage if T approaches T c . 
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FIG. 9. Incoherent c-axis conductance calculated for 
isotropic u = form factors, a momentum cutoff A = 7r/8, 
T = 0, J — 0.3, t' — —0.35 and three different dopings 8. 

We therefore have also studied finite values for A in the 
Gaussian in Eq.(7). Experimental evidence for strong 
forward scattering in the averaged squared tunneling 
matrix element has recently been found from the tem- 
perature dependence of the c-axis penetration depth in 



Y Ba2CuzO+ x [22]. Using the isotropic form factor u = 
and A = ir/8 Fig. 9 shows the conductance for three dif- 
ferent dopings at T=0. The spectra are dominated by a 
peak at approximately the frequencies 2V <& 2 + A 2 . This 
peak reflects the doping dependence of the "total" gap 
which increases monotonically with decreasing doping. 
In the upper panel the gap describes a superconducting 
gap, in the lower panel a d-CDW gap and in the middle 
panel a combination of both. The doping dependence of 
the main peak in Fig. 9 agrees well with experimental SIS 
spectra, see, for instance, Figs. 1 and 2 in [5], though 
the dips above the main line are more pronounced than 
in the experiment. Also negative conductances are only 
very rarely observed experimentally. In the slightly un- 
derdoped regime (middle panel in Fig. 9) the two order 
parameters coexist. The conductance shows in this case 
besides of the dominating high-frequency peaks a peak 
near the superconducting part of the gap. This peak 
moves towards smaller frequencies with increasing tem- 
perature and vanishes at T c . The exact energy position of 
this peak is somewhat below the superconducting part of 
the gap, 2A. This can easily be understood from Figs. 2 
and 3: Due to the interaction between the two gaps part 
of the BCS shoulder has been removed by the d-CDW so 
that the superconducting part of the gap appears smaller 
than the canonical value for 2 A. This reduced gap 
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FIG. 10. Incoherent c-axis conductance for five tempera- 
tures using a Gaussian distributed hopping matrix element 
with width A = tt/8. Increasing temperatures correspond to 
curves with energetically decreasing low-energy peaks. 

exhibits the expected temperature dependence as can be 
seen from Fig. 3. The upper and lower panels in Fig. 9 also 
show weak structures at low voltages. A closer inspec- 
tion, however, reveals that these structures are caused by 
the underlying band structure and are unrelated to the 
BCS gap. The spectra in Fig. 9 exhibit well-pronounced 
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dips at energies somewhat above the main peaks. These 
dips, which are also seen in tunneling spectra from break 
junctions, have been associated with self-energy effects 
due to the coupling to some boson [2,5]. We would like 
to stress that no self-energy effects have been taken into 
account in calculating Fig. 9. Simple model calculations 
indicate that one finds easily a dip in SIS spectra if the 
size of the region above 2A where spectral weight piles 
up because of the formation of the gap is comparable or 
smaller than the gap. 

Fig. 10 shows the temperature dependence of the inco- 
herent conductance for 5 — 0.114 and A = ir/8. The 
dominating higher peak is practically temperature inde- 
pendent for the temperatures shown in the figure. On 
the other hand, Fig.l indicates that both order param- 
eters vary in the considered temperature interval. One 
concludes from this that the higher peak reflects the to- 
tal gap which is rather independent of temperature. In 
contrast to that, the lower peak depends strongly on tem- 
perature. It moves towards zero frequency with increas- 
ing temperature, looses spectral weight and vanishes with 
vanishing A. Intrinis c-axis tunneling spectra in various 
cuprates show essentially the same features as in Fig. 10. 
In particular, the observed low-frequency peak, which 
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FIG. 11. Coherent c-axis conductance for anisotropic u=l 
form factors, T=0, J=0.3, t'=-0.35, rj = 0.01, and three dif- 
ferent dopings 6. 

seems not to be affected much by heating effects, also 
moves towards zero frequency with increasing tempera- 
ture and vanishes near T c . We would like to point out 
that this BCS-like peak can be seen in incoherent scat- 
tering only for a rather isotropic form factor 7. Other- 
wise, the tunneling of electrons near the nodal direction, 
where the SC gap is located, is too much suppressed. An- 
other prerequisite is that the averaged tunneling matrix 
element must be strongly momentum-dependent causing 



strong forward scattering. 

Besides of incoherent tunneling the quasiparticle cur- 
rent always contains a contribution J coh due to coherent 
tunneling, originating from the first term in the paran- 
theses in Eq.(5). J co h is given by Eq.(4) with g replaced 
by a(5k, q . The explicit expression for J co h thus becomes 

J coh (V) = ^^^(k) r Mm-f(" + eV)) 

,- J — 00 

k 

^2 ^1+21,1+2*' ( k, w + eV)A 1+2i 's+2i{K w). (13) 

ii'=0,l 

If the superconducting order parameter is zero, Eq.(l) re- 
duces to a 2x2 matrix. Calculating explicitly the spectral 
functions from this matrix and performing the frequency 
integration in the above integral one finds then that the 
sum over i, i' yields zero without any further approxima- 
tion. One thus obtains the important result that coherent 
tunneling is zero in a pure d-CDW. 

Fig. 11 shows the coherent conductance for u = 1, 
T = 0, and three different dopings. The spectral func- 
tions were obtained from the Green's functions using the 
frequency u + it] with 77 = 0.01. In the overdoped and 
slightly undcrdoped case the curves are similar to those 
in Fig. 9. In particular, corresponding curves have neg- 
ative conductances above the main peak caused by the 
restriction to small momentum transfers in the tunneling 
process. For S = 0.114 both curves show besides of the 
main peak associated with the total gap a second peak 
at smaller energies with BCS properties. However, one 
should note that Fig. 9 was calculated with u — whereas 
Fig. 11 with u = 1. It is somewhat surprising that coher- 
ent tunneling shows still the BCS peak though most elec- 
trons near the nodal direction are prevented from tunnel- 
ing due to the employed strongly anisotropic form factor. 
The dominance of small energy features in coherent tun- 
neling also is present in the underdoped case 5 = 0.077 
where the tiny BCS gap causes a sharp structure at very 
low energies. One important feature in Fig. 11 is related 
to the absolute values for the conductance depicted along 
the y-axis. The coherent conductance drops dramatically 
with decreasing doping. Going to smaller doping values 
one finds that J co h becomes zero if the superconducting 
order parameter vanishes in agreement with the above 
analytic result. Coherent tunneling is non-zero in case of 
a pure superconductor as shown by the upper panel in 
Fig. 11. However, it is zero for a pure d-CDW state and 
the d-CDW gap can only be probed in the presence of 
superconductivity. If the identification of the pseudogap 
phase with a d-CDW is correct coherent tunneling should 
vanish in the pseudogap phase. 

Fig. 12 illustrates the dependence of the coherent con- 
ductance on temperature for S = 0.114. With increas- 
ing temperature the position of the low-energy peak 
and its intensity decrease and approach zero at around 
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T = 0.0081 where the superconducting order parameter 
vanishes. The position of the high-energy peak as well 
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FIG. 12. Coherent c-axis conductance for anisotropic u=l 
form factors, J=0.3, t'=-0.35, rj = 0.01, and various temper- 
atures labeled according to decreasing maxima. 

as the dip are rather independent of temperature but the 
intensity of the whole high-energy part drops dramati- 
cally with temperature and vanishes at T — 0.0081. This 
again demonstrates that contributions from the d-CDW 
can only be seen in the conductance if the supercon- 
ducting order parameter is finite, i.e., in the coexistence 
regime. 

IV. SELF-ENERGY EFFECTS 

According to angle-resolved photoemission experi- 
ments the generic spectral function in the superconduct- 
ing state consists of a well pronounced peak followed by 
a dip and a hump towards larger energies [1]. In con- 
cordance with that the electron dispersion shows a kink 
between 30 and 70 meV below the Fermi energy [23]. 
These features occur throughout the underdoped, opti- 
mally doped and the overdoped regime. Most of these 
properties can be reproduced in a model where the elec- 
trons interact with a boson branch (which may be a 
phonon or a spin fluctuation) [24-26]. In the following 
we assume a dispersionless boson branch with a constant 
dimensionless coupling A. In the presence of SC and d 
CDW the inverse of the electronic 4x4 Green's function 
G(z, k) satisfies 

G- 1 (z,k)=G - 1 (z,k)-E(z,k), (14) 

where Gq 1 is given by Eq.(l) and £ is the self-energy. 
Because the boson-mediated interaction is momentum- 



independent in our model T, has only diagonal elements 
and it is En = £22 = £33 = £44, with 

_9^_ ( y> 2u s(E a )f(E a ) 

N ° t ^1 « z - ^) 2 - U^ a (E a - E P ) 
b(-u) )s(z - uj ) b{uj )s(z + u ) \ 

g 2 is related to A by g 2 = \iv / \2N (0)) , where N(0) is 
the density of states for one spin direction and loq the fre- 
quency of the boson. E a denote the four poles of G^°\ 
b the Bose function, and N c is two times the number of 
allowed momenta k. s(z) is given by 

s(z) = 2z(z 2 - (e 2 (k) + e 2 (k - Q))/2 

-A 2 (k)-$ 2 (k)). (16) 

Fig. 13 shows the real and imaginary parts of the re- 
tarded self-energy En(w + irj) at T = 0. The curves for 
different dopings look rather similar. The gap near the 
Fermi energy consists of the phonon energy plus the d- 
wave gap of the SC and/or d-CDW state. The panels 
in the figure illustrate the very smooth transition from a 
SC to a d-CDW gap with decreasing doping, passing also 
very smoothly through the coexistence region of SC and 
d-CDW. We used in this and in all the following figures 
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FIG. 13. Real (dotted line) and imaginary (solid line) parts 
of the self-energy for T=0, w = 0.065, 77 = 0.004, and A = 1. 

the value 1 for the dimensionless coupling constant A. 
This value corresponds to a change of the slope of the 
electron dispersion at the kink by a factor two in rough 
agreement with the photoemission data. 

Fig. 14 shows the SIN conductance using the u = 1 
anistropic hopping form factor and including self-energy 
effects. It is instructive to compare this figure with the 
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FIG. 14. SIN conductances with self-energy corrections for 
T=0, u=l, oj = 0.065, r? = 0.004, A = 1, and three dopings. 

analogous Fig. 7 where self-energy effects have been omit- 
ted. The BCS-structure seen in Fig. 7 for 5 = 0.114 has 
practically vanished in Fig. 14. The peaks in Fig. 14 are 
slightly broader, but the main effect of the self-energy is 
to move spectral weight from the main peaks to the side- 
bands. In the pure superconducting state at S — 0.178 
the sidebands consists of a clear dip and hump whereas 
in the two other cases the dip-hump feature is less pro- 
nounced. Both dip and hump move monotonically to- 
wards larger voltages with decreasing doping. The posi- 
tion of the dip in a pure superconductor is approximately 
half of the gap plus the boson energy. This rule also holds 
in the d-CDW and the mixed states. The exact differ- 
ences between the main peaks and the dip, however, fluc- 
tuate between 0.065 and 0.087 in Fig. 14. Though the dip 
in the SIN spectra is solely caused by the interaction with 
the bosons it may thus be difficult to determine precisely 
the boson energy from it. 

The solid and dashed lines in Fig. 15 are conductance 
curves for incoherent tunneling with and without self- 
energy effects, respectively. The bosons do not contribute 
to the non-diagonal self-energy because of the assumed 
momentum-independent coupling to the electrons. As a 
result, the bosons diminish both the SC and the d-CDW 
gaps via their diagonal self-energies. Consequently, the 
main peak moves towards lower frequencies but, consid- 
ered as a function of doping, this peak increases mono- 
tonically with decreasing doping as in the case without 
self-energy. Fig. 15 also illustrates that the distance be- 
tween the dip and the main peak is similar in the curves 
with and without self-energy effects and thus is rather un- 
related with the boson energy. For instance, at 5 — 0.178 
the distance between the dip and the main peak is 0.020 
and 0.018 in the case with and without self-energy, re- 



spectively, and thus much smaller than the boson energy 
of 0.065. 
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FIG. 15. Incoherent conductances with (solid line) and 
without (dashed line) self-energy corrections for T=0, u=l, 
A = oo, loq = 0.065, rj — 0.004, A = 1, and three dopings. 

If the experimental SIS spectra correspond to incoher- 
ent tunneling the boson energy cannot be obtained from 
the distance between the main peak and dip. Since the 
incoherent SIS spectrum is the folding of the SIN spec- 
trum in energy the incoherent SIS spectrum could have, 
in prinicple, a dip when the lower main peak is multi- 
plied by the dip at positive voltages and this dip posi- 
tion would be equal to the full gap plus the boson en- 
ergy. The upper panel of Fig. 15, however, shows that 
the folding in energy does not lead to a dip just below 
the maximum of the sideband. The main effect of the 
self-energy in Fig. 15 is to shift spectral weight from the 
main peak to the sideband consisting of a broad hump 
which monotonically moves towards larger voltages with 
decreasing doping. This hump is due to the folding of 
the (occupied) lower main peak with the (unoccupied) 
upper sidebands in Fig. 14. The solid and dashed lines in 
Fig. 16 represent coherent conductance curves with and 
without self-energy effects, respectively. We have omit- 
ted curves for the strongly underdoped case S = 0.077 
because they are smaller by two order of magnitudes due 
to the smallness of the supercondcuting order parameter 
in this case. Self-energy effects shift the main peaks to 
smaller energies, diminish somewhat the regions of neg- 
ative conductance, and create weak sidebands. Figs. 7 
and 14-16 suggest that self-energy effects and thus the 
nature of the boson spectrum appear more clear-cut in 
the SIN than in the SIS spectra. For instance, the dip 
in the SIN spectrum is soley caused by self-energy effects 
whereas that in the SIS spectrum is present even in the 
absence of any self-energy. The curves for the purely su- 
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perconducting case in Figs. 14 and 15 are similar to those 
published in Refs. [24,27]. 
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FIG. 16. Coherent conductances with (solid line) and 
without (dashed line) self-energy corrections for T=0, u=l, 
u>o = 0.065, f] — 0.004, A = 1, and two dopings. 



V. CONCLUSIONS 

The t — J model exhibits in the employed large-N limit 
a d-CDW phase at lower dopings besides of the supercon- 
ducting phase which is a natural candidate for the pseu- 
dogap phase observed in the cuprates. The density of 
states in the pure d-CDW state is strongly reduced near 
the Fermi level but still everywhere finite. This means 
that only part of the Fermi lines of the normal state are 
destroyed by the d-CDW and that the remaining Fermi 
lines form arcs around the nodal direction ending at the 
boundaries of the reduced Brillouin zone. With decreas- 
ing doping the length of the arcs become shorter. The 
ground state energy of the d-CDW can be lowered by 
introducing a d-wave superconducting gap near the arcs 
which explains the occurrence of a coexistence region of 
SC and d-CDW. Because the two gaps are well separated 
in k-space (the d-CDW gap resides near the antinodal, 
the superconducting gap near the nodal direction) fea- 
tures of the individual gaps survive even in the coexis- 
tence regime. 

In order to test the applicability of the above picture to 
cuprates we have calculated coherent and incoherent con- 
ductances and compared them with experimental spec- 
tra from break-junctions and intrinsic tunneling spec- 
troscopy. We find good evidence that the tunneling ma- 
trix element between layers is strongly anisotropic, sup- 
pressing tunneling of electrons near the nodal direction, 



which is in agreement with band structure arguments. 
Incoherent tunneling thus probes mainly electrons near 
the maximal gap at the X and Y points. This gap trans- 
forms in a very smooth way from a superconducting gap 
at large dopings to a d-CDW gap at small dopings pass- 
ing continuously through the coexistence regime. Cal- 
culated incoherent conductances thus fit best to the ob- 
served spectra from break-junctions which are character- 
ized by one peak moving monotonically to larger voltages 
with decreasing doping. We find that coherent tunnel- 
ing is only non-zero for a non-vanishing superconduct- 
ing order parameter. Below T c it shows in spite of the 
anisotropic tunneling matrix element two peaks which 
can be associated with SC and d-CDW. The appearance 
of a low-energy peak in the calculated coherent tunneling 
spectrum, which moves to if T c is approached from be- 
low, is unique for intrinsic tunneling spectra. From this 
we conclude that tunneling in stacked, intrinsic junctions 
is dominated by coherent tunneling and that the appear- 
ance of the low-energy peak related to superconductivity 
supports models with two competing order parameters 
in the underdoped region. Including self-energy effects 
due to the coupling of electrons to a dispersionless boson 
branch as suggested by ARPES removes part of the re- 
gions of negative resistances and also creates sidebands 
which, at least in the case of SIN junctions, resemble 
those which have been measured. 

The authors thank Secyt and the BMBF ( Project 
ARC 99/007) for financial support and A. Yurgens and 
V.M. Krasnov for useful discussions. 
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